State transition of a non-Ohmic damping system 
in a corrugated plane 
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Anomalous transport of a particle subjected to non-Ohmic damping of the power 5 in a tilted 
periodic potential is investigated via Monte Carlo simulation of generalized Langevin equation. It 
is found that the system exhibits two relative motion modes: the locking state and the running 
state. Under the surrounding of sub-Ohmic damping (0 < 8 < 1), the particle should transfer into a 
running state from a locking state only when local minima of the potential vanish; hence the particle 
occurs a synchronization oscillation in its mean displacement and mean square displacement (MSD). 
In particular, the two motion modes are allowed to coexist in the case of super-Ohmic damping (1 < 
5 < 2) for moderate driving forces, namely, where exists double centers in the velocity distribution. 
This induces the particle having faster diffusion, i.e., its MSD reads (Ax 2 (t)} — 2D^t SeB . Our 
result shows that the effective power index S B s can be enhanced and is a nonmonotonic function of 
the temperature and the driving force. The mixture effect of the two motion modes also leads to a 
breakdown of hysteresis loop of the mobility. 

PACS numbers: 05.70.Fh, 05.60.-k, 05.40.-a, 05.10.Ln 



I. INTRODUCTION 

There are many physical situations which can be de- 
scribed by Brownian transport in a tilted periodic poten- 
tial, for example, Josephson junction [l[, charge-density 
wave 0, superionic conductor rotation of dipoles in 
external field [1] , phase- locking loop [|| , diffusion on sur- 
face [f| , and separation of particles by electrophoresis Q • 
The quantitative properties of those systems have been 
discussed in plenty of literatures [8| , such as the depen- 
dence of the coherence level of transport on the tem- 
perature, driving force and shape of potential @; the 
huge enhancement of the effective diffusion coefficient 
relative to the free diffusion [l(| ITl|. Jl3| : the response 
of output to the noise and signal [ljand so on [l4l |. 
Since theoretical tools and numerical algorithms are not 
sufficient in non-Markovian dynamics with a frequency- 
dependent non-Ohmic damping, most of the models are 
established in the Ohmic damping environment. How- 
ever, the frequency-dependent damping is more general 
because a large number of stochastic processes fail to be 
the Markovian dynamics. 

Recent studies on anomalous diffusion and transport 
are mostly limited in the absence of potential, linear force 
case or the sub-diffusion in a potential [l5[ . It is worth to 
point out that the behavior of a particle moving in a pe- 
riodic potential immersed in the super-Ohmic damping 
environment might be far more complicated than that in 
the Ohmic and sub-Ohmic damping cases. In compar- 
ison with the previous findings for great enhancement 
of the effective diffusion coefficient floL [ill [l2j and the 
hysteresis loop of mobility [8| in the Ohmic damping en- 
vironment, we will perform a detailed investigation in the 



present work on the diffusion and the mobility of a parti- 
cle subjected to arbitrary non-Ohmic damping in a cor- 
rugate plane. This is in terms of an effective algorithm 
proposed by us [l6| to numerically solve a generalized 
Langevin equation (GLE) with arbitrary damping kernel 
function and corresponding thermal colored noise. 

The paper is organized as follows. In Sec. II, we de- 
scribe briefly the anomalous transport model by means 
of the GLE. In Sec. Ill, the two basic quantities of in- 
terest: the generalized effective diffusion coefficient and 
the fractional mobility are defined; the novel behaviors 
of diffusion and mobility are shown. Finally, we draw a 
conclusion of our findings in the section IV. 



II. THE MODEL 

We consider a Brownian particle moving in a one- 
dimensional periodic potential under the influence of non- 
Ohmic memory friction and a constant external driving 
force. The dynamics of the particle is governed by the 
following GLE [13, El, 

x(t) = v(t), (1) 
mv(t) = -ml j(t - t')v(t')dt' + U'(x) + y/mk B T£(t), 
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where ks is the Boltzmann constant, T is the temper- 
ature of the environment, j(t) is the damping kernel 
function and related to £(t) through the fluctuation- 
dissipation theorem [13, El] 

mt(t'))=7(\t-t'\), (2) 

where is a zero mean Gaussian colored noise and its 
spectral density reads 

a^lvUMy-VM). ,3, 



2 



The small |w| behavior of (|£(w) 2 |) 1S a power-law char- 
acterized by the index (5 — 1. The function f c (\u)\/u) r 
a high frequency cutoff function of typical width oj c 
and u> -C ui c denotes a reference frequency allowing for 
the constant mjs to have a dimension of viscosity for any 
5 [H| • The cases of < 5 < 1 and 1 < S < 2 are the sub- 
Ohmic damping and super-Ohmic damping, respectively; 
S = 1 is the Ohmic one, then 7(0;) is equal to a constant 
and the noise is white. In Eq. (1), U(x) is considered to 
be a tilted periodic potential, 



? (.,;) = — f7 cos ( 



Fx. 



(4) 



The minima of U (x) vanish when the driving force F is 
taken to be the critical value: F c = Uq2ti/\ = 1.0. 

In the calculation, the natural unit (to = 1 and fc^ = 
1), the dimensionless parameters: Uq = 1.0, A = 2ir, 
js = 4.0, the smooth cutoff function f c = exp(— lo/uj c ) 
[2l| with uj c = 4.0, and the time step At = 0.01 are used. 
The test particles start from the origin of coordinate and 
have zero velocity, here 2 x 10 4 test particles are used to 
describe the stochastic distribution of a Brownian parti- 
cle. 



III. DIFFUSION AND MOBILITY 

The quantities of foremost interest are the diffusion 
coefficient and the mobility. Here we generalize the both 
into non-Ohmic damping case with an arbitrary power 
index 5, 



D (S) := lim 



1 



t^oo T(l + <5) 



D 5 t (Ax 2 (t)) s , 



V* : = , lim Trr-7s^oD s t (x{t)) s , 



(5) 
(6) 



where $D t denotes the fractional derivative. The algo- 
rithm for numerically calculating the two quantities is 
presented in the Appendix [X] 



A. Diffusion 

We begin our studies from the situation of sub-Ohmic 
damping. We have noticed that the case of sub-diffusion 
dynamics has been discussed by Goychuk and Hanggi 
[22], where the GLE and the fractional Fokker-Planck 
equation approaches to the escape dynamics are used 
and compared. The escape is governed asymptotically 
by a power law whose exponent depends exponentially 
on both, the barrier height and the temperature. If the 
ratio of the barrier height to temperature is too large, 
the diffusion motion in a washboard potential well below 
a critical tilt cannot be observed numerically within a 
reasonable time window, i.e., nearly all the test particles 
are confined in the locking state. Therefore, we consider 
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FIG. 1: (a) Time-dependent MSD of the sub-Ohmic damp- 
ing particle of 5 = 0.6 and (b) its mean displacement. The 
parameters used are: F = 5.0 and T — 0.1. 



the transport of sub-Ohmic damping particle subjected 
to a large driving force because of the efficiency of nu- 
merical simulation of GLE. Figures 1 (a) and 1 (b) show 
time-dependent MSD and MD of the particle of S — 0.6. 
With increasing the driving force F until local minima 
of the potential vanish, we find that the MSD of the par- 
ticle shows a quasi-periodic oscillation. As long as the 
MSD of the particle experiences a quasi-period process, 
the particle will move the distance of a periodic length 
Ax = 2tt along the direction of external force. 

In Fig. 2, we plot the space probability distribution 
of a sub-Ohmic damping particle at different times. It 
is seen that with the evolution of time, the width of 
the probability distribution becomes periodically narrow 
when the particle moves in the bottom of a potential 
well; the one is broad when the particle arrives at the 
top of potential. Unlike the normal Ohmic damping par- 
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FIG. 3: The MSD of the particle with S = 1.7 at a large 
driving force F — 5.0 for various temperatures. 



tide [lOj, its probability distribution is not centralized. 
Our results can be interpreted as follows. Under the sub- 
Ohmic damping environment, the particle has a strong 
memory to its initial position and thus the diffusion in 
the coordinate space is much slow. If the potential have 
local minima, it is quite difficult for the particle to escape 
from the potential well, thus the particle is in a locking 
state during the period of simulation. As the local min- 
ima of the potential vanish, the particle subjected to a 
large driving force can enter the running state. There- 
fore, its distribution width is modulated by the periodic 
structure of potential and thus the MSD of particle oc- 
curs a quasi-periodic oscillation. 



Figure 3 shows that the quasi-periodic oscillation phe- 
nomenon becomes unconspicuous when the temperature 
increases for T > 1.0 at F = 5.0. In this case the parti- 
cle transfers completely into the running state and hence 
the structure of the potential might have less influence 
on the transport process. 

In particular, in the case of super-Ohmic damping, we 
find numerically that the asymptotic MSD of the particle 
can be approximately written as a power function, 

(Ax 2 (t))=2D^t s ^ T ' F \ (7) 

where S e g depends on T and F. Indeed, the index 5 c s is 
not always equal to S like the Ohmic and the sub-Ohmic 
damping cases, but varies non-monotonically with F . For 
a moderate F, we find a prominent result: The effective 
power index 5 e s exceeds 2 (i.e., the ballistic diffusion [23|, 
EH, [HI) when the periodic potential is titled observably 
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FIG. 5: The space probability distribution of the particle at 
time t = 50.0. The inset figure is the probability distribution 
at the locking state only. The parameters used are: S = 1.7, 
F = 0.75 and T = 0.1. 



but its local minima still exist. Further analysis shows 
that the mysterious diffusion behavior is caused by the 
mixing of the locking state and the running state. 

In Figs. 4 (a) and 4 (b), we plot the effective power 
index 5 e g as a function of F for various T and 6. It is 
seen from Fig. 4 (a) that the maximal value of 6 c s versus 
F forward with increasing temperature. Similar behavior 
of d e ff for other super-Ohmic damping cases can also be 
observed, as shown in Fig. 4 (b). The smaller the value 
of S is, the larger F where the maximum (5 e ff appears. 
This can be interpreted qualitatively as follows. In the 
non-Markovian rate theory [2^] , for a sufficiently big ra- 
tio of barrier height to temperature, the super-diffusion 
in a periodic potential should turn out into the normal 
diffusion because the escape events are exponentially dis- 
tributed in time and no overlong jumps can occur, so that 
6 c fi — 1 when the tilt of the potential is small. Also, for 
very large F and the structure effect of potential vanish- 
ing, all the test particles can move into a running state, 
thus (S c ff = 5. However, for a middle tilt, some test par- 
ticles are confined in the potential well (in the locking 
state) and others drift quickly forward (in the running 
state). A proper proportion between the locking state 
and the running state should induce the maximum of the 
effective power index. Apparently it becomes easier for 
the particle in the locking state to escape the well and 
join the running state with the increase of either T or S. 
Therefore, the maximum of S c g appears in the case of a 
small F at high temperature; in the case of large S at low 
temperature, respectively. 

In Fig. 5, we illustrate the case of S = 1.7 at the low 
temperature (T = 0.1) and the middle tilt (F = 0.75) 
to depict the coexistence of the two motion modes. The 




FIG. 6: (Color online) The velocity distributions of the par- 
ticle at t = 150 (a) and t = 330 (b). The solid and 
dashed lines are the results of the supper-Ohmic S — 1.7 and 
Ohmic 5 = 1.0 cases, respectively. The parameters used are: 
F = 0.75 and T = 0.1. 



backward and forward barrier heights are given by 
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For a low temperature T < U2 <C U\, the particle oscil- 
lates around the potential well, however, it still escapes 
over a barrier with a small probability. Since the bar- 
rier crossing process is quite slow at low temperature, the 
particle in the locking state has an approximate Gaussian 
space distribution centering at xq [xq — ^ arcsin( 2 ^ )], 
as shown in the inset of Fig. 5. Once the test particle 
climbs over the barrier, it will no longer be restricted 
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again. Because the next hill of the titled periodic po- 
tential is lower than the present one, the kinetic energy 
of the particle gaining from the external driving force is 
greater than the dissipated energy due to the memory 
friction. The particle can enter the running state, so it 
drifts quickly along the direction of driving force. This 
results in a long tail appearing in the space probability 
distribution of the particle. For the Ohmic damping par- 
ticle, once escaping over a barrier, it will slide to the next 
well and be trapped again. Hence the space probability 
distribution of the Ohmic damping particle is smashed 
into some small pieces and does not make up a running 
state. 

Figures 6 (a) and 6 (b) show the coexistence of the 
two motion modes of the super-Ohmic damping particle, 
which can be more intuitively for the velocity distribu- 
tions of the particle at times t = 150 and t — 330. As 
expected, we do not hnd the coexistence phenomenon of 
two velocity modes appearing in the normal diffusion. It 
is seen that the super-Ohmic damping particle with an 
increasing probability enters the running state; the differ- 
ence between the two center velocities in the running and 
rocking states becomes large with the evolution of time. 
This implies that as long as the test particles escape out 
of the well, they will be accelerated by the corrugated 
plane and join the running state. Of course, the present 
locking state occurs simply because we cannot observe 
the motion on the time-scale of our numerical simulation 
as the escape rate are very rare [22], but a new locking 
state relative to the forward running state should arise 
once the original locking state disappears. 



B. Mobility 
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The mobility determined by Eq. (6) as a function of 
the driving force is plotted in Fig. 7. Starting from zero 
tilt and switch adiabatically on the tilt F, the mobility 
of the particle remains zero when all the test particles 
are in the locking state until some of them join the run- 
ning state. In comparison with the Ohmic damping case, 
we find that the hysteresis loop is broken and become 
stagger. At that point A, the mobility jumps to infin- 
ity and then drops to a constant at point B and keeps 
this constant with increasing F (point C). The point B 
corresponds to the critical force F = F c , where the lo- 
cal minima of the corrugated plane vanish. When the 
driving force decreases adiabatically, all the test parti- 
cles are kept in the running state and the mobility ap- 
proaches a constant until the driving force becomes so 
small that most of the test particles are trapped in the 
potential wells. At the point D, the mobility falls to zero 
again. While for the usual Ohmic damping case shown 
in Fig. 7 (b), we find that a bistability occurs in the 
region 0.7 < F < 1.4, which forms a hysteresis loop. 
However, for the sub-Ohmic damping case, we have not 
found similar hysteresis loop of mobility as what arises 
in the Ohmic damping case. 



FIG. 7: (Color online) (a) The fractional mobility times the 
damping constant as a function of driving force in the super- 
Ohmic damping case with 8 — 1.7. (b) The normal Ohmic 
damping result with 7 = 0.5. The solid and dashed lines cor- 
respond to the forward and backward processes, respectively. 
The temperature is T = 0.1. 



IV. SUMMARY 

We have investigated the transport of a non-Ohmic 
damping particle in a titled periodic potential and re- 
ported a prominent finding: The titled periodic poten- 
tial as a simple equipment which not only enhances the 
diffusion coefficient, but also changes diffusive behavior 
of the particle. This is due to a novel phenomenon of 
two motion modes: the locking state and the running 
state, which can appear and transfer in the corrugated 
plane. In the sub-Ohmic damping case, the mean square 
displacement of the particle shows a quasi-periodic prop- 
erty when the driving force is larger than the critical 
value where the minima of the potential vanish. While for 
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the super-Ohmic damping case, the two motion modes, 
namely, there exists two centers in the velocity distribu- 
tion, can be coexisted and transferred. Thus the power 
index for the mean square displacement of the particle is 
enhanced. In comparison with the hysteresis loop of mo- 
bility of the normal case, the hysteresis loop of mobility 
of a super-Ohmic damping particle is broken. 

The anomalous Brownian motion in a periodic po- 
tential is representable for many applications occurring 
in areas such as in condensed matter physics, chemical 
physics, molecular biology, communication theory, and so 
on. We are confident that both theoretical and experi- 
mental works in the future will help one to clarify further 
and shed more light onto all these intriguing issues and 
problems. 
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APPENDIX A: NUMERICAL METHODS FOR 
FRACTIONAL CALCULUS 

The so-called Riemann-Liouvill fractional integral is 
defined through [27[ 

t Ii 5 f(t) = ^ £ dt> - [W_ s , t > 0, S > 0, (Al) 

whereas its left-inverse t Df reads 

t Df := D™t I t s - m ,m- 1< 6< m,m € N, (A2) 

where oD™ denotes the ordinary derivative of order m. 
In this present work, the case of to = is concerned. For 
completeness, we define 

J t ° = D° = I, (A3) 

where I is the identity operator. It is convenient to make 
use of the discrete operators of translation (shift) and fi- 
nite differences to derive the approximative recursive ex- 
pressions to the fractional differentiation operator oDf. 
In a lucid way the theory of numerical differentiation and 
integration (with equidistant grid points) has been devel- 
oped in C hap ters 7 to 10 of Ref. [28|. See also Chapter 
6 of Ref. 

Let r £ E, we define the shifting operator E T and the 
backward difference operator V T by them acting on a 
function u(t) for <€l, 

E T u(i) = u{t + r), (A4) 

V T u(t) = u(t) - u(t - t). 
We furthermore have the relation, with I as the identity 
operator, 

V T =I-£T r . (A5) 



Using these notations, we write the approximation [u(i) — 
u(t—r)]/h for the derivative u'(t) of a differentiable func- 
tion u(t) as Vh,u(t)/h for a small positive h with accuracy 
0(h) as the function u(t) is sufficiently smooth. High or- 
der derivatives u^ nS> (t) =o D™u(i) ( n € N) with small 
h > 0, can be approximated by 

[W^u(t)]/h n = h- n (I - E->Yu(t), (A6) 

again in case of u(t) being sufficiently smooth, with or- 
der of accuracy 0(h). The powers Vi n ' can be readily 
expanded via the binomial theorem 

3=0 VJ ' 

This leads to the known formula 

h- n ^(-\y(^\u(t-jh) = D?u(t) + O(h). (A8) 
j=o ^ 3 ' 

The remarkable fact now is that these formulas can be 
generalized to the case of non-integer order of derivative. 
Replacing the positive integer n by a positive real number 
6 amounts to use the formal power 

V£ = f>l)^W^, (A9) 
3=0 ^ J ' 

in analogy to the expansion (E~ h replaced by the com- 
plex variable z) 




(A10) 



which is convergent if \z\ < 1. We thus get the Gfunwald- 
Letnikov approximation: 

h- 5 V{u(t) = h- s Y,(-iy( ■ Mt-jh) = D s t u(t)+0(h). 
3=0 ^ J ' 

(AH) 

If u(t) decays to zero sufficiently fast as t — > oo , in 
particular if u(t) = for t < 0, will not diverge and 
hence for the latter case, we have 

[t/h] 

h- s V s h u(t) = h- s J2 (-I)''! «(* - 3h). (A12) 

3=0 ^ J ' 

By using the property of Gamma function 

mm -5) = -^— 

sm(7rd) 

we obtain the final recursion of calculating fractional cal- 
culus: 

oDfu(kh) = ^ £ ^r—YM(k j)h). (A13) 
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